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A MULTILOOP IMPROVEMENT OF NON-SINGLET QCD EVOLUTION EQUATIONS 

S. V. Mikhailov |] 

Joint Institute for Nuclear Research, Bogoliuhov Laboratory of Theoretical Physics, 

14 1980, Moscow Region, Dubna, Russia 

An approach is elaborated for calculation of "all loop" contributions to the non-singlet evolu- 
tion kernels from the diagrams with renormalon chain insertions. Closed expressions are obtained 
for sums of contributions to kernels P{z) for the DGLAP equation and V{x, y) for the "non- 
forward" ER-BL equation from these diagrams that dominate for a large value of 60, the first 
/3- function coefficient. Calculations are performed in the covariant ^-gauge in a MS-like scheme. 
It is established that a special choice of the gauge parameter ^ = —3 generalizes the standard 
"naive nonabelianization" approximation. The solutions are obtained to the ER-BL evolution 
equation (taken at the "all loop" improved kernel), which are in form similar to one- loop solu- 
tions. A consequence for QCD descriptions of hard processes and the benefits and incompleteness 
of the approach are briefiy discussed. 

PACS: 12.38.Cy, 12.38.-t, 13.60.Hb 



I. INTRODUCTION 



^ Evolution kernels are main ingredients of the well-known evolution equations for the parton distribution of DIS 
S^ • processes [|l| and for parton wave functions yj in hard exclusive reactions. These equations describe the dependence 

jrt of parton distribution functions and parton wave functions on the renormalization parameter ^ . The calculations 
performed beyond the one-loop approximation for the forward DGLAP evolution kernel P{z) P,H, and what is 
more, for the nonforward Efremov-Radyushkin-Brodsky-Lepage (ER-BL) kernel V{x^y) [gJ6) were challenged and 
complicated technical tasks. 15 years later, the 3-loop results for these kernels are not known yet, except for the 
first few elements of anomalous dimension in DIS, obtained numerically in [^]. In this situation, it seems useful 
to try other ways to gain knowledge about high-order corrections to these kernels and to the solutions to the 
corresponding equations. 

Here I discuss the results of the diagrammatic analysis and multiloop calculations of the DGLAP kernel P{z) 
and ER-BL kernel V{x^y) in a certain class of the "all-order" approximation of perturbative QCD (pQCD). The 
corresponding diagrams include the chains of one-loop self-energy parts (renormalon chains) into the one-loop 
diagrams (see Fig. 1). The regular method of calculation and resummation of the indicated classes of diagrams 
for these kernels based upon their simple forest structure has been suggested in B. There was established that 
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the resulting series possesses a nonzero convergent radius, therefore the infrared rcnornialons are absent in the 
kernels. The results of that summation for both the kinds of kernels (DGLAP and ER-BL) obtained earlier in the 

Nf 

framework of a scalar model in six dimensions with the Lagrangian Lint — 9 /, (VA*V'i¥')(6) '^ith Nf of the scalar 

i 

"quark" flavours (jjji) and "gluon" (ip) are analyzed here for non-singlet QCD kernels. For the readers convenience 
some important results of the paper g would be recalled. 

The insertion of the chain into the "gluon" line ("chain-1" in 0) of the diagram in Fig.l a,b and resummation 
over all bubbles transforms the one-loop kernel aPo{z) — az = a(l — z) into the "improved" kernel P^^'{z; A) 



aPo{z)^az''^~'P^'\z;A) = az 



^-A,, .^7^(0)^ 



(z)--(l-A)- 



where A = aNf^^{0), a = . ,3 . (1) 



7^ (A), 

Here, 7^(ip)(e) are one-loop coefficients of the anomalous dimensions of quark (gluon at Nf — 1) fields in D- 
dimension (D = 6 — 2e) discussed in g; for the scalar model 7^(e) = J,p{s) = B{2 — e, 2 — e)C(e), and C(e) is a 
scheme-dependent factor (C(0) = 1) corresponding to a certain choice of an MS-like scheme. The argument A of the 
function 7i^(A) in (|l|) is a standard anomalous dimension (AD) of a "gluon" field. On the other hand, the result (|l|) 



(-^)"n„r.i , o/o^" „„.„., „ .2 (-^) 



n-2 



Nn-2 



corresponds to resummation of a class of seriea^ like a ; — (lii[^] + 8/3)"-series, aA — — — (Infzl)" -series, 

'-' n! [^ ^ 2)! 

. . ., (see Table 1 in g]) into the kernel which dominate at large Nf. 

The resummation of this "chain-1" subseries into an analytic function in A should not be taken by surprise. 

Really, the considered problem can be connected with the calculation of large Nf asymptotics of ADs' in order of 

1/A^/. An approach was suggested by A. Vasil'ev and colleagues at the beginning of the 80'es g to calculate the 

renormalization- group functions in this limit, they used the conformal properties of the theory at the critical point 

9 — dc corresponding to the non-trivial zero gc of the D-dimensional /3- function. This approach was extended by J. 

Gracey for the calculation of ADs' of composite operators of DIS in QCD in any order n of PT, jl^]. I used another 

approach which is close to [O; contrary to the large Nf asymptotic method, it does not appeal to the value of 

parameters NfTji, Ca/'2 or Cp, associated with different kinds of loops in QCD. To illustrate this feature, let us 

consider the insertions of chains of one- loop self-energy parts into the "quark" line of diagram Fig. la ("chain-2" 

in g). Contributions of these diagrams calculated in the framework of the above scalar model do not contain the 

parameter Nf, nevertheless, they can be summarized into the kernel p(^^(z; B) S 



aPoiz) = az ^"^-' p(^^\z;B) =az{l + B^ 



,,.-B 7v-(0) 



; where B = 07^(0), z = 1 — z, (2) 



according to the same approach. This corresponds to summation of various series like (n+l)-a ; — Infz] -\ — 

. . . into the kernel. The operator (1 + _B d/dB) appearing in front of formula (y) expresses an inherent combinatoric 
factor (71 + 1) to these diagrams. Following that line, the "improved" QCD kernel p(^^(z; A) was obtained in iQ 
for the general case of a mixed chain (quark and gluon bubble chain) in ^- gauge. 

Here, we present the QCD results similar to Eq.(|l|), in the covariant ^- gauge for the DGLAP non-singlet kernel 
P{z;A). Analytic properties of the function P{z;A) in variable A are analyzed. The assumption of the "Naive 
Nonabelianization" (NNA) approximation ||l^ for the kernel calculation ||lj| is discussed and its generalization 
based on ^ = —3 gauge is suggested. The numerical importance of the resummation in this case is demonstrated. 



^On the other hand these class simply corresponds to the Taylor expansion of the kernel P''^\z; A) in a new parameter A, 
so, the n-term of expansion corresponds to the n-bubble chain insertion. 



The ER-BL evolution kernel V{x,y) is obtained in the same multiloop approximation as the DGLAP kernel, by 
using exact relations between the P and V kernels [g,g] for a class of "triangular diagrams" . The considered class 
of diagrams represents the leading 6o-contributions to both the kinds of kernels. Partial solutions for the ER-BL 
equation, <i>„(a;,y4), are derived. The multiloop "improved" kernels P{z;A)^ V{x,y;A) and solutions ^n{x,A) are 
compared with the exact results in 3(2)-loop approximation. 

II. TRIANGULAR DIAGRAMS FOR THE DGLAP EVOLUTION KERNEL 



Here, the results of the bubble chain resummation for QCD diagrams in Fig.l for the DGLAP kernel are discussed. 
These classes of diagrams generate, in particularly, the contributions ^ as {A ln[l/z])" /n\ in any order n of pQCD. 
Based on the resummation method of Ref. |8| in the QCD version, one can derive the kernels ^(1"^'^^^) (corresponding 
to the diagrams in Fig.l a,b,c) in the covariant ^— gaugCQ whose explicit expressions are presented in |12|| . They 
contribute to the total kernel p(^^(z; A, ^) that has the expected "plus form" 



P^^\z;A,0^asCF2- 
asPoiz) = asCp^ ■ 



zz-''{l-AY + 



2z 



A{A,0' 



(3) 

(4) 
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FIG. 1. The diagrams in figs, la - Ic are "triangular" diagrams for the QCD DGLAP kernel; dashed lines 

for gluons, solid lines for quarks; black circles denote the sum of all kinds of the one-loop insertions (dashed 

circles), both quark and gluon (ghost) or mixed chains; the slash on the line denotes the delta function S{z — kn) 

{k is the momentum on the line) which is traced to the representation of the composite operator ®, see [6] for 

details; MC denotes the mirror-conjugate diagram; Id is an example of a diagram for the nonforward ER-BL 

kernel. 

a 1 

where as = —^, Cp — {N^ — l)/2Nc, Ca = Nc, T^ — - are the Casimirs of SU(A^c) group, quantity 
47r 2 

A = A(0, ^) = — as7g(0, ^), and the 7g(0,^) is the one-loop coefficient of the standard AD of the gluon field. For 
comparison with (j^, the one-loop result asPQ{z) is written also down, the latter can be obtained as the limit 
P^-'^^(z;A -^ 0,^). The function yl(e,^) is defined as A{e,^) = — as7g(e,^), where the function 7g(e,C) is the 



'The gauge parameter ^ is defined via the gluon propagator in the lowest order iD^,j{k ) = 



fc2 -|- ie 



9i^v 



(«-l)^ 



coefficient of the anomalous dimension in D-dimcnsion, here and below D = 4 — 2e. In other words, it is the 
coefficient Zi{e) of a simple pole in the expansion of the gluon field rcnormalization constant Z that includes both 
its finite part and all the powers of the e-expansion. So, one can conclude that the "all-order" result in (^ is 
completely determined by the expression 7g (e) (73 (e,C)) ^'^^ the single quark (or/and gluon) bubble subgraph. 
The function 7(,(£,0 thus defined is an analytic function in the variable e by construction, see ||]. Note that the 
function 0^73 (e,^) at zero reveals itself again as the argument of the same function, Eq.(P), through the variable 
A. The quantity A here plays the role of a new perturbative expansion parameter. Equation (0) is valid for any 
kind of insertions, i.e., jg — 7^ for the quark loop, -fg = -fg^ for the gluon (ghost) loop, or for their sum 

l,iAO=li'\A)+^l^\A,0; (5) 

when both the kinds of insertions are taken into account. The 6{l — z) - terms appearing in the partial contributions 
(see [|l2[) are exactly accumulated in the form of the [...]+ prescription in (||), and the ^ - terms successfully cancel. 
This is due to the evident current conservation for the case of quark bubble insertions; including the gluon bubbles 
into consideration merely modifies the effective AD 'jg' (A) -^ 7g(A, ^), conserving the structure of result (^, see 
PJl^. Substituting the well-known expressions of ^g{e) from the quark or gluon (ghost) loops (see, e.g., |l5[ ) 

^l^He) = -SNfTnBiD/2,D/2)Cie), (6) 

7^g''\s,0^^BiD/2-l,D/2-l) 




(7) 

into formula (0) one can obtain P^^'{z; A, ^) for both the quark and gluon loop insertions simultaneously. Here, the 
coefficient C(e) = r(l — £)r(l-|-e) implies a certain choice of the MS scheme where every loop integral is multiplied 
by the scheme factor T{D/2 — l){fj,'^ /AnY (MSi scheme). The rcnormalization scheme dependence of P^^'>{z;A) 
is accumulated by the factor C(e). For another popular definition of a minimal scheme, when a scheme factor is 
chosen as exp(c ■ e), c = —je + • ■ • instead of T{D/2 — 1) (MS2 scheme), the coefficient C(e) does not contain any 
scheme "traces" in final expressions for the renormalization-group functions. 

Of course, the final result (|3|) will be gauge-dependent in virtue of the evident gauge dependence of the gluon 
loop contribution 7g (£,0- ^ '^"w expansion parameter A in this case. 



A^ 



sls{0,O = -as{i'H0,O + l^fH0) 



''^'-'^^C.~UfTn 



(8) 



^3 2 J 3 

is the contribution to the one-loop rcnormalization of the gluon field. The positions of zeros of the function jg{A, ^) 
in A, which represent the poles of P{z;A,^), also depend on S,. The kernel P'-^-'(z; A,^) becomes gauge-invariant 
if we restrict themselves only to the quark-loop insertions, i.e., jg — > 7^ , A — > A^"^' — —asj^'^'{0) = ag—TjiNf, 
and P^^^(z; A, ^) is reduced to Pq '{z;A^'^^), as it is presented in |S]. It is instructive to outline analytic properties 
of Pq '{z;A^'i^) in A^"^^ based on Eq. (^) and on the explicit form for jg' in (|^): (i) the range of convergence of 
the PT series corresponds to the left zero of the function 'yg'{A) and is equal to Ao = 5/2, which corresponds 
to a^ = 15n/Nf, so, this range looks very broad ||, Ug < bn at Nj — 3; (ii) the resummation into Pq {z;A) is 



^Here we consider the evolution kernel P{z,A) itself. We do not consider that the factorization scale fj? of hard processes 
would be chosen large enough, fj? > m^, where the p-meson mass nip represents the characteristic hadronic scale. For this 
reason, the used coupling as{fJ.'^) could not be too large. 



substantial, two zeros of the function Pq {z; A) in A appear within the range of convergence (in MSi scheme). 
Of course, the moments of this reduced kernel Pq (z; A^'^^) agree with the generating function for the anomalous 
dimensions, obtained earlier in [[lO|. 



III. A MODIFIED NNA VERSION FOR KERNEL CALCULATIONS 

The expansion of Pq (z; A) in A provides the leading as ■ (cLgNj ln[l/z])" dependence of the kernels with a large 
number Nf in any order n of PT [p|. But these contributions do not numerically dominate for real numbers of 
flavours Nf —4, 5, 6. That can be verified by comparing the total numerical results for 2- and 3-loop ADs' of 
composite operators (ADCO) presented in with their A^/-leading terms, see Table 1. There the contributions to 
coefficients of different Casimirs in the ADCO are presented. To obtain a satisfactory agreement at least with the 
two-loop results, one should take into account the contribution from next-to-leading A^/-terms. As a first step, 
let us consider the contribution from the completed renormalization of the giuon line, which should generate a part 
of next-to-leading terms. Below, we examine an exceptional choice of the gauge parameter ^ = —3. For this 
gauge the coefficient of one-loop gluon AD 7g(0, —3) coincides with bo, the one-loop coefficient of the /3-function Qj 
and A = —asbo. Therefore this gauge can be used for reformulating the so-called |1^] NNA proposition to kernel 
calculations. Note, just this value of S, has been used in M& to estimate the total gluon contribution only from 
the gluon bubble in order a^ to the process of e+ e~ annihilation. Other interesting applications of this gauge to 
approximate the exact loop results have been considered in |17| , p^ . 

To obtain the NNA result in a usual way, one should substitute the coefficient bo for 7g (0) in the expression 
for a'*) by hand (see, e.g., uM). Note, the use of such an NNA procedure does not improve Pq {z;A) and leads 
to poor results even for the two-loop level, i.e., for the aj. Pi{z) term of the expansion, see M]. The NNA trick 
expresses common hope that the main logarithmic contribution can follow from the renormalization of the coupling 
constant gs. The first effective realization of this idea goes back to the well-known Brodsky-Lepage-Mackenzie 
(BLM) prescription for the scale setting [Q formulated in the next-to-leading order approximation. That gs- 
renormalization appears as a sum of contributions from all the sources of renormalization of gs at the vertices of 
triangular diagrams. Let us consider the gluonic, vertices, and quark line renormalizations successively in the case 
of the ^ = —3 gauge. The one-loop gluon renormalization in this gauge imitates the contributions from all other 
sources and the coefficient bo appears naturally via of 7c,(0, — 3). At the same time, in the one-loop vertices 
renormalization constant Zip, 



1-ZiF 



C^e + %(3 + 



the nonabelian part vanishes at ^ = —3, while the corresponding Abelian part, UgCp^, is compensated by the 
renormalization of the quark line of a triangular diagram, —QsCf^, due to the abelian Ward identityfl. So, due 
to the cancellations, only the gluon contribution survives in gs renormalization and provides the expected 6o-terni, 
as6olii[2;]. These properties of cancellation can be illustrated by the well-known diagram by diagram results for 
two-loop Pi{z) presented in Feynman gauge in ^^ (for Vi{x,y) in ||j2l[]). Indeed, the terms, connected with the 



Here, for the /?(a3)-function we adapt f3{ag) = —boa^ — bi a^ . . ., bo = — Ca Nf, bi = — C^ — Nf I 2Cf H Ca I , • 



quark field/ vertic rcnormalization are proportional to ln[l — z] in these diagrams and really cancel in the gauge 
invariant sum of all contributions. In contrast to that, the ln[z]-ternis collect the coefficient asbo. Though we 
should not take into account the self-energy chain ("chain-2" in the Intr.) and "rainbow" graph insertions into the 
quark line unless the vertices of the triangular diagram, dressed in the same manner, is included into consideration, 
we see that their contributions should be cancelled in the first log-parts for the discussed gauge. For these reasons 
we can guess the gauge ^ = —3 "exceptional" for the one-loop chain dressing. 

To analyze the resulting effect of "all-loop" resummation for the case ^ = —3 in M), let us choose the common 
factor A/A{A, —3) in formula (0) (below the notation a ~ agbo = —A is introduced). 



P'^^\z--a,~3)^asCF2- 



zz"(l-Fa)2 + 



2z 



l+a 



A 



&o 



r(2 + 2a)(3-t-2a) 



1-z 
bo 



A 



A{-a,-3)' 



(9) 
(10) 



A{~a,-3) jgi-a,-3) (r(l + a))^ C(-a) (4CAa2 + a(36o + 2Ca) + 35o) ' 

for a crude measure of the modification of the kernel in comparison with the one-loop result asPo(z). The factor 
(as well as the whole kernel P{z; —a, —3)) has no singularity in a for a > 0. Considering the curve of this factor in 
the argument a in Fig. 2, one can conclude: 
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FIG. 2. The curves of the factor A/A{~a, —3) in a; the solid line corresponds to MSi scheme; the dashed line corresponds 
to MS 2 scheme 



(i) the factor 



A 



noticeably grows with argument a in the range of the standard PT validity. Really, 



A{-a,-3) 

this factor reaches 1.32 for the MSi scheme (1.17 for the MS2 scheme), if we take the naive boundary of validity 

An 
of the standard PT, oq = 0.5, a^o — — ao ~ 0.7 that corresponds to the value of as on the hadronic scale; thus, 

°° . . 
the resummation is numerically important in this range, see Fig. 2. 

(ii) scheme dependence looks not too strong for acceptable values of parameter a. 

Note that Eqs.(||, ^j could not provide the valid asymptotic behavior of the kernels for z ^ 0. A similar z-bchavior 
is determined by the double- logarithmic corrections which are most singular at zero, like as (os In [z]) |22| ]. These 
contributions appear due to rcnormalization of the composite operator in the diagrams by ladder graphs, etc.. 



rather than by the triangular ones. But, Eq.fl) can provide a main z-behavior for not too smaU z due to simple- 
logarithmic corrections. To obtain the low boundary of this z-region, let us compare effects from simple and double 
logarithmic contributions taking into account the main singular terms up to 3 loops; 



Uio 



P{z) = as2CF 



-al2CF [■■■ + bo \n[z] + (2Ca - 3Cf) In^H] 



a%2CF 



llfeo - 2 , M^^o, 2. A , CI-%\Cf-CaI2\\ 4. 
, ^ ln[z] + — In {z\\ + ^ In [ 



(11) 



The first terms in the square brackets in ( pd| ) follow directly from the expansiorjw| of Eq.(^ in a; the second term 
in the second line is the double-log from the exact two-loop calculations; and the last term in the third line was 
predicted by J.Blumlein&A.Vogt in p2]. From (O) rough estimate follows to the boundary of validity of Eq.(Q), 



z ~ 0.1 — 0.05 at moderate as — 0.3 — 0.1. The most singular In [z]-tcrm in (O) becomes important for z < 10"'^. 
It seems naturally to combine the improved by the simple-logs kernel P^^'{z; —a, —3) with the first double-logs 
contribution from the exact two- loop calculations pj into a modernized kernel P{z), 

P{z) = P^^\z- -a, -3) + alCp [{Po{z)Ca - (1 + z)Cf) \n [z] - 4{Cf - Ca/2)Po{-z)F{z)] ^ , 

which works up to z ~ 10"'^. 

At the end let us consider the integral characteristics of the kernel P'-^^z; —a, —3) to compare with the exact 
results. The expansion of this kernel in a generates partial kernels aj.Pti\{z), a^Pt2)iz)y ■ ■ ■ which in turn produce 
ADCO a^ r(i)(n), a^ r(2)(7^), . . • according to the relation r(n) = /„ dzz"P{z). Let us compare these elements of 
ADCO and a few numerical exact results from m collected in Table 1: 

(i) evidently, the leading A^/-contributions are reproduced exactly for any T/j\{n); 

(ii) we consider there the next-to-leading A^/-contributions to the coefficient r(i)(n) generated by gluon loops and 
associated with the Casimirs CfCa/2, the C|n~term is missed, but its contribution is numerically insignificant. It 
is seen that in this order the C^C^-terms are rather close to exact values (the accuracy is about 10% for n > 2) 
and our approximation works rather well; 

(iii) in the next order, the contributions to r(2)(n.) associated with the coefficients Nf ■ CfCa and C\Cf arise, 
while the terms with the Casimirs C%, Nf ■ Cp, C^Ca are missed. The involved Casimirs [Nf ■ CfCa, C\Cf), 
as it follows from the estimations in j^], also dominate numerically in the third order ADCO r(2)(n), which gives 
a hint for success. Nevertheless, contrary to the previous item, all the generated terms are opposite in sign to 
the exact values, and the "^ — —3 approximation" doesn't work at all. So, we need the next step to improve the 
agreement with 3-loop results - to obtain the next-to-leading iV/-terms by an exact calculationQj. 



'"'"The expansion of the Eqs.lMBI) in the next orders generates the Rieman zeta-functions started with ("(3) in order of 0{a*) 
§§Here the double-Zog's part is rewritten from @; F{z) = \ \n^[z\ - 21n[z] ln[l + z\- 2Li2{-z) - Li2{l) 
***An example of similax calculation in QCD recently has been demonstrated in p3] 



Table 1. The results of Fj-i 2) (") calculations ( r(n) = L dxx"P{z)) performed in different ways, exact numerical 
results from H and approximation obtained from P{z,A,£^) with ^ — —3; both numerical and analytic exact 
results are marked by the bold print. 





r(i)H 


r(2)H 


CfCa 


NfCp 


C^Cf 


Nf ■ CfCa 


Nf ■ Cf 


n=2 












Exact 


13.9 


-2.3704 


86.1 + 21.3 C(3) 


-12.9- 21.3 C(3) 


-0.9218 


e = -3 


11.3 




-42.0 


12.9 




n=4: 












Exact 


23.9 


-4.9152 


140.0+ 19.2 C(3) 


-18.1 -41.9 C(3) 


-1.5814 


e = -3 


23.5 




-76.0 


23. 




n=6 












Exact 


29.7 


-6.4719 


173 + 19.01 C(3) 


-20.4- 54.0 C(3) 


-1.9279 


e = -3 


31.1 




-95.6 


28.5 




n=8 












Exact 


33.9 


-7.6094 


196.9+ 18.98 C(3) 


-21.9- 62.7 C(3) 


-2.1619 


e = -3 


36.3 




-109.0 


32.3 




n=10 












Exact 


37.27 


-8.5095 


216.0+ 18.96 C(3) 


-23.2 - 69.6 C(3) 


-2.3366 


e = -3 


41.00 




-119.28 


35.24 




n=12 












Exact 


40.02 


-9.2555 


7 


7 


-2.4753 


e = -3 


44.64 




-127.61 


37.58 





On the other hand, it looks rather naive to expect a good agreement of the values of r(2)(n), obtained from 
expansion of (H) with exact three-loop results. Really, our approach takes account only of the simplest two-bubble 
chain diagrams (to be precise, only 3 types of diagrams) among all the set of 3-loop diagrams. What can we 
expect from the next-to-leading Nf corrections, like a^ ■ (as(asiV/)"), mentioned above in item (iii)? The set of 
corresponding diagrams starts with a part of all 3-loop diagrams, and this part, as I hope, dominates at this loop 
level. Let us consider the diagrams underlying the "tower" of these corrections; these diagrams contain: 



(a) only one- loop insertion into giuon lines or vertices, the number of such diagrams in a covariant gauge amounts 
to 39 (without M.C. diagrams); these diagrams can be obtained from the set of two-loop diagrams presented, e.g., 
ini; 

(b) essentially two-loop self-energy insertions into gluon lines, they are only of 3 types of the diagrams presented 
in Fig. l(a,b,c); now the black bubble denotes the sum of these two- loop self-energy parts. 

The calculation of the contributions from diagrams of type (a) looks as a formidable task; the result is substantial 
and could not be guessed a priori. On the contrary, diagrams of type (b) lead to partially an expected contribution 
to the kernel, 

asPo{z) ■ al^'^^\OHz) + ■ ■ ■ , (12) 

which is evidently generated by the two-loop AD, a^^g (^), of the gluon line in Fig. 1, where 



{0 = -rCi-NA2Cp + -CA\-\^\ (e-l)U+^ ' (13) 



(see, e.g., |24[| ). It seems tempting to include that contribution into consideration via a modification of basis 
formula (|3|) , even though "by hand" . Namely, substituting a new "corrected" expansion parameter A* for A 

A^A* = -a,7,(0,0-a'7^''(6 

into expression (|^), one can restore the contribution ( [l2[ ) in the expansion of this model P'^^^z; A* , £) kernel. As 
the next step, one should choose a new "corrected" value of the gauge parameter, ^* = —3 + 0{ag) -^ ^. Following 
the NNA idea and our previous reasoning about an exceptional ^ = —3 gauge, let us define it by a natural condition 
through the /^-function 

7s(0,r) + a,l\^\n = ^'o + ash, + 0{al), (14) 

that leads to the value 

r = -3 + as^-CA [nj - ^-C^ + 0{al). 

The hypothesis on P^'^\z; A* , (^*) only slightly reduces the discrepancy between the exact and model results for 
r(2)('^) in Table 1. Moreover, it generates the contributions to a new required Casimir Nf ■ Cp, which appear of 
the same sign, and are in order smaller than the exact ones. It is clear, that the model P'^^^(z; A*,^*) is a step 
along the right direction, but it is obviously insufficient. So, we insist on accurate calculations of both the types 
(a) and (b) diagrams to obtain a reasonable approximation to exact r(2)(r7,)-results. 

IV. THE NONFORWARD ER-BL EVOLUTION EQUATION AND ITS 

SOLUTION 

Here we present the results of the bubble resummation for the ER-BL kernel V{x, y). The latter can be derived in 
the same manner as it was done for the DGLAP kernel P{z), see Appendix A in M. On the other hand, V{x, y) can 
be obtained as a "by-product" of the previous results for P{z), i.e., we use again [|8|,p^ the exact relations between 
the V and P kernels established in any order of PT Q for triangular diagrams. These relations were obtained 
by comparing counterparts for the same triangular diagrams considered in "forward". Fig. la, and "nonforward" , 
Fig. Id, kinematics. 



Collecting the contributions from triangular diagrams, see [n2J, one arrives at the final expression for V^-'-'^ in 
the "main bubbles" approximation 



y«(x,y;A0=asCF2 



\ 1-^4 / 1 



e{y>x)[-] [I- A 



yj \ y-x 



ff..^.„-., ,., 



that has a "plus form" again due to the vector current conservation. The contribution V'^'^^ in (|l5|) should dominate 
for Nf ;2> 1 in the kernel V . Besides, the function V^^^^^ {x, y; A, ^) possesses an important symmetry of its arguments 
X and y. Indeed, the function V{x,y; A,£^) = V^^'{x,y; A,£^) ■ {yy)^~ is symmetric under the change x <-> j/, 
V(x, y) = V{y, x). This symmetry allows us to obtain the eigenfunctions ijjn{x) of the "reduced" evolution equation 

ill 

1 

V'^'\x, y- A)My: A)dy = r(n; A)V„(x; A), (16) 



V'„(2/; A) = [yyY^'^^^--^ N{nA) ' ^^'^^ ^^^"^^ = ^^^ ~ ^^Z^' ^^ = * - 2^' (1^) 



N{n,A)=2^-^''^^''\T{n + 2d^{A))/(n\ [n + d^{A)) {V{d^{A)f)) 



y{a)/ 



and d^ip{A) is the effective dimension of the quark field when the AD A is taken into account; C„ (z) are the 
Gegenbauer polynomials of an order of a; N{n, A) is the norm of C," {y — y), |25|. The partial solutions <i>(a;; Og, /) 
of the original ER-BL-equation ( where I = ln(/i^//iQ)) 

d 



A^ ..,2 



2-^$(a;;a„0= / V^''){x,y; A) <P{y;as,l)dy (18) 



are proportional to these eigenfunctions tpn{x;A) for a special case of the stopped evolution as — aj, /3(aJ) = 0, 
see, e.g., [£6|,||. The result ( [17| ) for the eigenfunctions at ^ = —3, has been confirmed in |g^ by "a partial 
resummation of conformal anomalies" and in a suggestion of a large value of bo. Let us examine ipn{x', —a) in ( |l7| ) 
as an approximation to the exact two-loop solution derived in a closed form in p6|. Expanding, e.g., "ipoiy, ~a) in 
parameter a we can express ipQ^'''^{x) versus the exact solution ?/;g^'"^*(x) 

-00 (a;; -a) -^ ^^'''"'(x) = 6xx < 1 + 0^60 ( ln(a;x) + o ) f ' (1^) 

0^^"^*(:e) =6xxil + asbo (hxx) + ^^ + a^Cp (a^ (^) + 2 - y) | • (20) 

The term 'ijJq^^^ [x) coincides with the "conformal symmetry-predicted" (CSP) part in (EQ), ( proportional to bo), 



this part dominates in ?/;g^°'^* (x) in the mid-region of the parameter x, 0.3 < a; < 0.7. The other part in ( |2C| ) is 
generated by the "additional conformal symmetry breaking term" pq] ; it contributes in the opposite phase to the 
first one and it is large and enhanced near the end points. For the latter reason, ip'^PP^{x) become useless at n > 2 
even for the mid-region x description, see [ p6[ . 

In the general case P{as) ^ let us start with an ansatz for the partial solution of Eq.([lq), ^nix;as,l) ^ 
Xnio-sil) ■ i^n{x;A), with the boundary condition Xri(as,0) — 1; <I>„(x;as,0) ~ ^pn{x;A). For this ansatz. Eg. ([l8|) 
reduces to 

{ti^d,,2 + /3(a,)aaJ In {^n{x; as, I)) = r(n; A). (21) 

In the case n = 0, the AD of the vector current r(0; A) — 0, and the solution of the homogeneous equation in ( pl| ) 
provides the "asymptotic wave function" 



$o(x; as, I) = M^; A) - ^^(^-^((1 - x)x)^'-^^ (22) 

where A — — as(/i^)7(0, ^), and a.sifJ.'^) is the running coupUng corresponding to a /3-function P{as). Similar solutions 
have been discussed in [ |l4| in the framework of the standard NNA approximation. Solving simultaneously Eq. ( pi| ) 
and the renormalization-group equation for the coupling constant a^, we arrive at the partial solution $„(x;as,^) 
in the form 

{j-a-sifj.'^) Yin A) 1 
- / / da } . (23) 

Ja,i^.l) /3(a) J 

An adequate choice of /3-function in ( p3| ) must correspond to the same modified NNA approximation that was 
applied for r(n, A) calculation, but it is absent yet. The /3-function in a large iV^ expansion, that is equivalent to 
quark bubbles resummation, has been computed in p8||. 



V. CONCLUSION 

Here, I present closed expressions in the "all order" approximation for the DGLAP kernel P{z) and ER-BL kernel 
V{x, y) resulting from resummation of a certain class of QCD diagrams with the renormalon chain insertions. The 
contributions from these diagrams, P'^^'(z; A) and V'^^\z;A), give the leading Nf dependence of the kernels for 
a large number of flavours Nf ^ 1. These multiloop "improved" kernels are generating functions to obtain 
contributions to partial kernels like aV P(„)(z) in any order n of the perturbation expansion. Here yl ~ a^ is 
the new expansion parameter that coincides (in magnitude) with the anomalous dimension of the gluon field. On 
the other hand, the method of calculation suggested in [p| does not depend on the nature of self-energy insertions 
and does not appeal to the value of parameters NfTfj, Ca/2 or Cp associated with different loops. This allows 
us to obtain contributions from chains with different kinds of self-energy insertions, both quark and gluon (ghost) 
loops, see [n2|. The price for this generalization is the gauge dependence of final results for P'-^^z; A{^) , S,) and 
V^'^''{z; A{S,), £,) on the gauge parameter ^. 

The result for the DGLAP nonsinglet kernel P^^\z;A,^) is presented in (||) in the covariant ^-gauge, it looks 
similar in form to the simple one-loop kernel. The analytic properties of this kernel in the variable as are dis- 
cussed for an exceptional gauge parameter ^ = —3. This choice of the gauge allows one to generalize the naive 
nonabelianization suggestion and provides the leading feo-behavior of the kernel for large bo ^ 1. For this gauge 
P^^^{z; A, —3) in (^) works up to z ~ 0.1 — 0.05 at moderate as — 0.3 — 0.1, and reproduces two-loop anomalous 
dimensions a^r(i)(n) with a good accuracy, while the standard "naive nonabelianization" proposition fails at this 
level. But on the next three loop level the "^ — —3 approximation" is insufficient, see quantities Vi2\(n) in Table 
1. At the end, a hypothesis about a possibility to extend the approach to 3-loop level is briefly discussed. 

The contribution V'^^^ix, y; A, ^) to the nonforward ER-BL kernel (nS) is obtained for the same classes of diagrams 
as a "byproduct" of the previous technique jlj,^ . The partial solutions ( p^ ) , ( [2^ ) to the multiloop improved ER-BL 
equation are derived, that are similar in form to the one-loop solutions. The form of these solutions appearing at 
f ~ —3 was confirmed independently in [p7| . The lowest harmonic ipoi^i ^) roughly imitates the cc-behavior in the 
mid-region of the exact two-loop solution ( [p6[). 

The obtained results are certainly useful for an independent check of complicated computer calculations in higher 
orders of perturbation theory, similar to ||7|] ; they are useful for the analysis of evolution "at small x" ; they may be 
a starting point for further multiloop approximation procedures. 
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